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Abstract. We study the reversible quantum Stirling cycle with a single spin or two coupled spins as the 
working substance. With the single spin as the working substance, we find that under certain conditions 
the reversed cycle of a heat engine is NOT a refrigerator, this feature holds true for a Stirling heat engine 
with an ion trapped in a shallow potential as its working substance. The efficiency of quantum Stirling 
heat engine can be higher than the efficiency of the Carnot engine, but the performance coefficient of 
the quantum Stirling refrigerator is always lower than its classical counterpart. With two coupled spins 
as the working substance, we find that a heat engine can turn to a refrigerator due to the increasing 
of the coupling constant, this can be explained by the properties of the isothermal line in the magnetic 
field-entropy plane. 


1 Introduction 

In thermodynamics, a cycle consists of a collection of ther¬ 
modynamic processes transferring heat and work, while 
changing pressure, temperature, and other state variables, 
eventually returning a system to its initial state. In the 
process of going through this cycle, the system may per¬ 
form the conversion of heat or thermal energy to mechan¬ 
ical work, named as a heat engine. Heat engines distin¬ 
guish themselves from other types of engines by the fact 
that their efficiency is fundamentally limited by Carnot’s 
theorem. There are four basic thermodynamical processes 
in classical thermodynamics, i.e., isothermal process, iso- 
choric process, adiabatic process and isobaric process, 
they correspond to fixed temperature, volume (the gen¬ 
eralized coordinate of the system), entropy, and pressure 
(the generalized force conjugated to the generalized coor¬ 
dinate), respectively. For example, an isothermal process 
is a change of the system, in which the temperature re¬ 
mains constant. This typically occurs when the system 
is in contact with an thermal reservoir, and the change 
occurs slowly enough to allow the system to continually 
adjust to the temperature of the reservoir through heat ex¬ 
change. In contrast, an adiabatic process is where the sys¬ 
tem exchanges no heat with its surroundings. The working 
substance can be any system with a non-zero heat capac¬ 
ity, then one may wonder: what happens with a quantum 
system as the working substance. These questions were an- 
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swered in references [1-6], in which four classical processes 
have been generalized to the quantum case. With these 
generalizations, all the classical thermodynamical cycles, 
such as Carnot cycle, Otto cycle, Brayton cycle, Diesel 
cycle, and Stirling cycle [7], can be generalized to the 
quantum case. The working substance for a quantum heat 
engine or quantum refrigerator can be all kinds of quan¬ 
tum systems, ranging from single spin [8-11], two-level 
or multilevel systems [12-16], harmonic-oscillators [17] to 
coupled spins [18-20]. Based on these working substance, 
the properties for quantum Carnot cycle [21-23], quan¬ 
tum Otto cycle [19,24], quantum Brayton cycle [20,25], 
and quantum Diesel cycle [6] have be considered from 
different points of view including the positive work con¬ 
dition [3,4,15], efficiency and power [26-28]. Apart from 
the pure academic interest, the quantum generalization of 
those thermodynamical processes may find applications. 
In fact, several schemes have been proposed to realize the 
quantum heat engine, some of them are realized in exper¬ 
iments. For example, driven coupled spins linked to two 
split heat baths can serve as a quantum Carnot heat en¬ 
gine by adjusting the energy difference [29]. The quantum 
Carnot cycle may also be realized by a single spin cou¬ 
pled to an oscillator and sandwiched between two heat 
baths by designing the dissipative process [30]. Recently, 
experiment to realize quantum Otto cycle are presented by 
confining a single ion in a tapered linear Paul trap and by 
coupling the ion to engineered laser reservoirs [31], or by 
cold bosonic atoms confined to a double well potential [32] . 
The main motivations about this topic can be summa¬ 
rized as follows. First, with the development of quantum 
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information processing and nanotechnology, one may ma¬ 
nipulate and control a system at the nanoscale. At the 
nanoscale, quantum effects dominate. This invokes people 
to consider the thermodynamics of small quantum sys¬ 
tems. Second, the study of the quantum heat engine pro¬ 
vides us a way to better the work extraction and a method 
to improve the efficiency of heat engine beyond the clas¬ 
sical limit [33,34]. Quantum resources, such as quantum 
coherence and entanglement [35], are resources which have 
no classical counterpart. It has been reported that the use 
of these quantum resources can strikingly affect the ther¬ 
modynamics. For example, we can extract work from a 
single heat bath by vanishing quantum coherence [36-38], 
we can obtain a heat engine with efficiency higher than 
the Carnot one [39] by virtue of squeezed heat bath. Last 
but may be the most important, these studies provide us 
with a new viewpoint to better our understanding of the 
fundamental concept, such as Maxwell’s demon and the 
universality of the second law [40-42], and the third law 
of thermodynamics [43,44], in thermodynamics, statistical 
physics and quantum mechanics [45-52]. 

In the previous works, the focus is on the cycles where 
there are two quantum adiabatic processes. There exists 
cycles in classical thermodynamics that have no adiabatic 
process. Then a question arises: what will happen when 
this type of cycle is generalized to the quantum domain. In 
this Manuscript, we study the quantum Stirling heat en¬ 
gine and refrigerator, which is an extension of the Stirling 
cycle from classical thermodynamics to quantum physics. 
An ideal Stirling cycle consists of two isochoric processes 
and two isothermal processes. The Stirling cycle is re¬ 
versible, which means that if supplied with mechanical 
power, it can function as a heat pump for heating or cool¬ 
ing. We will show that there are many interesting proper¬ 
ties which are different from its classical counterpart. For 
example, with a single spin-| as the working substance, 
we show that under certain conditions, the reverse cycle 
of a heat engine is not a refrigerator, because it absorbs 
heat from the bath at high temperature and releases heat 
to the bath at low temperature, while the work is done 
on the working substance. This feature can be observed 
in the heat engines with an ion trapped in a shallow po¬ 
tential as its working substance. In addition, we find that 
with two coupled spins as the working substance, a heat 
engine can be turned into a refrigerator by increasing the 
coupling constant. 

2 Single spin 

A quantum Stirling cycle is a quantum counterpart of clas¬ 
sical Stirling cycle, which consists of two quantum isother¬ 
mal processes and two quantum isochoric processes. In 
the quantum isothermal process, the working substance is 
kept in contact with a heat bath at a constant tempera¬ 
ture T and the system remains in thermal equilibrium with 
the heat bath at every instant of time. The generalized co¬ 
ordinate of the working substance is changed slowly, as a 
result, the energy spectrum E n and the occupation prob¬ 
abilities p n , as well as the internal energy U are changed. 



Fig. 1 . Scheme illustration of a quantum Stirling heat engine 
cycle based on magnetic field-energy B-U diagram. 


For a quantum system at thermal equilibrium, the occu¬ 
pation probabilities p n are given by p n = exp(— (3E n ), 
where Z is the partition function, H is the Hamiltonian 
of the system, f3 = with k the Boltzmann constant. It 
should be noted that in the classical isothermal process, 
where the ideal gas is taken as the working substance, 
the internal energy does not change, but it is not true for 
quantum isothermal process. In classical isochoric process, 
the volume V, i.e., the generalized coordinate is fixed. Ac¬ 
cordingly, in a quantum isochoric process, the generalized 
coordinate of the working substance is fixed. As a result, 
the energy spectrum of the system keeps unchanged dur¬ 
ing the quantum isochoric process and no work is done. 
For the detailed discussions about these generalizations, 
we refer the reader to reference [5]. 

The Hamiltonian of the working substance under con¬ 
sideration can be written as H s = \Bg z , where o z is the 
usual Pauli operator and B is the external magnetic field. 
There is only one parameter B in the Hamiltonian (i.e., 
the generalized coordinate), hence the quantum isochoric 
process for the working substance means that the mag¬ 
netic field keeps as a constant and the magnetic field ex¬ 
erting on the working substance is a quantum analogy of 
the volume in the classical working substance (e.g., ideal 
gas). With these knowledge, the four stages for the heat 
engine can be depicted as follows. Stage 1 : A — > B is an 
isothermal process. In this process, the system interacts 
with an heat bath at high temperature TV The magnetic 
field changes from B 2 to B\ slowly such that the system is 
at thermal equilibrium with the heat bath at any instance 
of time. Heat is absorbed from the bath (named Q 1 ). 
Stage 2: B —» C is a quantum isochoric process. In this 
process, the magnetic field keeps B\ and the temperature 
of the system decreases from Th to T c . No work is done 
and heat is released (named Qbc) hr this process. Stage 3: 
C — » D is another isothermal process. The magnetic field 
changes from B\ to B 2 . The temperature of the system 
keeps T c . Heat is released to the heat bath (named Q 2 ). 
Stage 4 • D —> A is another isochoric process. Heat is ab¬ 
sorbed (named Qda) and there is no work. A schematic 
diagram for this heat engine is given in Figure 1. All the 
above heat exchanges can be calculated as follows. In the 
two isothermal processes, we have dQ = TdA, where S 
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is the entropy of the working substance. For a two-level 
system at thermal equilibrium with temperature T, the 
entropy is: 
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In the two quantum isochoric processes, no work is done. 
The heat transfer equals to the change of the internal 
energy U. For our system, the internal energy reads: 
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As a result, Qb c and Qda are: 

Q bc = U(C) - U(B) = B (tanh tanh ^° Bl 


Qda = U(A) — U ( D ) = -A (tanh T 2 — tanh ^ h 2 


( 5 ) 


( 6 ) 


Fig. 2. Heat engine efficiency rj as functions of B\ (in units 
of kT c ) and B 2 /B 1 . Other parameters are kTh = 3 kT c . Here 
rj = | is the classical Carnot efficiency under this condition. 


Here we have three possible cases: (i) AQ = 0. This is 
the condition for perfect regeneration. As discussed above, 
this is true for classical ideal gas as the working substance. 
With quantum mechanical systems as the working sub¬ 
stance, however, a restriction on the two magnetic fields 
Hi, F> 2 , Th and T c is necessary, namely, these parameters 
can not be chosen arbitrarily. Usually the cycle falls into 
the following two categories ((ii) and (iii)). (ii) AQ < 0. In 
this case, the amount of heat Qb c flowing from the working 
substance to the regenerator in one regenerative process 
is larger than that (i.e., Qda) flowing from the regener¬ 
ator to the working substance in the other regenerative 
process. This requires that the redundant heat in the re¬ 
generator must be released timely to the cold heat bath, 
otherwise the regenerator may not operate normally. The 
release of heat for the regenerator increases the amount 
of heat release to the heat bath T c . (iii) AQ > 0. In con¬ 
trary to the case (ii), the amount of heat flowing from the 
working substance to the regenerator Qb c is smaller than 
that Qda flowing from the regenerator to the working sub¬ 
stance. The inadequate heat in the regenerator must be 
compensated from the heat bath Th- This results in the 
increasing of heat absorbed from heat bath at high tem¬ 
perature. Hence the heat absorbed from heat bath Th per 
cycle can be obtained as: 


where Q > 0 and Q < 0 correspond to absorption and re¬ 
lease of heat, respectively. Based on these heat exchanges, 
we can obtain the work output per cycle according to the 
first law of thermodynamics as: 


W = Y,Q = kT h In 
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In classical Stirling heat engine, a regenerator is used in or¬ 
der to improve the performance of the cycle in the two iso¬ 
choric processes. The regenerator is applied to the quan¬ 
tum Stirling heat engine too. This means that the heat 
released in Stage 2 is sent back in Stage 4. As a result, 
Stirling cycle is usually categorized as regenerative cycle. 
For classical ideal gas, we always have \Qb c \ = \Qda\ (or 
Qbc + Qda = 0). However, it is not the case for the cycle 
with quantum working substance. Hence a net amount of 
heat transfer between the system and regenerator in the 
two isochoric processes is needed as: 



Qh — Q 1 + 5AQ, (9) 


where 5 = 0 for AQ < 0 and S = 1 for AQ > 0. 

We first discuss the positive work condition for quan¬ 
tum Stirling cycle. It has been shown that [3,4] the positive 
work condition of quantum Otto heat engine with a sin¬ 
gle spin as the working substance is Th > |^-T C , which is 
more restrictive than that in the classical Otto heat en¬ 
gine. Here we can prove that the positive work condition of 
quantum Stirling heat engine with a spin-^ particle as the 
working substance is Th > T c and B\ < £> 2 , which is con¬ 
sistent with the classical one. The proof for this condition 
is given in Appendix A. 

Next we consider the heat engine efficiency. It can be 
obtained as: 

W 



The analytical expression of the engine efficiency is com¬ 
plicated, to show it clearly, we plot the numerical result in 
Figure 2. For classical Stirling heat engine, the efficiency is 
rj = 1 — j?-, which is the same as that of Carnot efficiency 
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r]c and is independent of the generalized coordinate of the 
system (volume for P — V — T system). With quantum sys¬ 
tem as the working substance, this observation does not 
hold. We observe from the figure that, the efficiency de¬ 
pends on both B\ and B 2 . In the figure, we use the ratio of 
B 2 to B\ (namely, B 2 /B 1 ) as the axis. This ratio can be 
seen as an isothermal compression in quantum systems. 
In most regime, there is an optimal isothermal compres¬ 
sion for a fixed B\, the efficiency reaches the peak, which 
corresponds to the condition for perfect regeneration. In 
other words, when AQ = 0, the efficiency arrives in its 
maximum and this maximum is usually higher than the 
Carnot efficiency r]c . We should note that this result does 
not violate the second law of thermodynamics, as such a 
high efficiency profit is obtained from the effects of the 
regenerator. In order to work normally, the regenerator 
needs to consume additional energy. If we remove the re¬ 
generator, the heat engine cycle can still operate normally. 
However, the efficiency would decrease and it is lower than 
the Carnot efficiency. We should emphasize that the effi¬ 
ciency for a classical Stirling cycle with optimal regenera¬ 
tor equals to the Carnot efficiency. Our results show that 
the use of the regenerator can improve the efficiency of a 
quantum thermodynamical cycle higher than the Carnot 
efficiency. 

In the following we focus on the inverse cycle of the 
Stirling heat engine. We should note that our heat engine 
is reversible. The proof is given in Appendix B. In the 
frame of classical thermodynamics, it is proved that the 
inverse cycle for a heat engine is a refrigeration cycle. How¬ 
ever, this is not the case for quantum heat engine. Similar 
to the heat engine cycle, we start from point A in Fig¬ 
ure 1. A —> D is an isochoric process. The heat exchange 
Qad = — Qda • Other processes can be described similarly 
with Q c 5 = — Qbc , Q'i = —Q 1 , Q f 2 = —Q 2 * The schematic 
diagram is the same as in Figure 1 but all arrows point 
to opposite direction as in Figure 1. The total work can 
be obtained as W' = —W. According to the positive work 
condition discussed above, we have W' < 0 when T c < Th 
and B\ < B 2 . This indicates that the work is done on the 
working substance. The net amount of heat transfer be¬ 
tween the system and the regenerator in the two isochoric 
processes is AQ' = Qad+Qcb • When AQ' < 0, the amount 
of heat Q a d flowing into the regenerator in one regenera¬ 
tive process is larger than that ( Q c b ) flowing out from the 
regenerator in the other regenerative process. The redun¬ 
dant heat in the regenerator must be released to the heat 
bath T c . This results in the reduction of the net amount 
of heat absorption from the hot bath T c . Hence the total 
heat absorption from cold bath per cycle is: 

Qc = Q' 2 -S\AQ'\ ( 11 ) 


where S = 1 for AQ' < 0 and S = 0 for AQ' > 0. We can 
see that if AQ' < 0 and Q' 2 < \AQ'\, we have Q c < 0. This 
means that the heat is released to the cold bath. In this 
case, the inverse cycle is not a refrigerator. In the case of 
Qc ^ 0 ? the coefficient of performance for the refrigerator 
is 



£ 



0 1 2 7-> 3 4 5 


Fig. 3. Coefficient of performance £ as functions of B± (in 
units of kT c ) and B 2 /B 1 . Other parameters are kTh = 3 kT c . 


Figure 3 shows the dependence of the coefficient of perfor¬ 
mance on B\ and the isothermal compression B 2 /B 1 . In 
this figure, the blank region corresponds to W' < 0 and 
Qc ^ 0* As we discussed, the cycle is not a refrigerator in 
this regime. The condition that the cycle in this regime 
is Q' 2 + AQ' < 0. In addition, the coefficient of perfor¬ 
mance depends on B\ and £> 2 , which is also different from 
the classical result £ = T . Compared with Figure 2, 
we find the maximum stands for the condition for perfect 
regeneration, i.e., when AQ' = 0, the coefficient of perfor¬ 
mance reaches its maximum and this maximum is usually 
lower than the classical counterpart. 

To shed more light on the condition with that the in¬ 
verse cycle of a Stirling heat engine is not a refrigeration 
cycle, we consider the case of f3h = 0 and B\ — 0, which 
means that the temperature of the hot bath is very high, 
and the energy of the lower level of the spin is zero. With 
these conditions, we have: 


Q c = — kT c In cosh((3 c B 2 /2) < 0, (13) 


which is equal to zero when (3 C B2 = 0. This is illustrated 
by the blank region very close to the horizontal axis in 
Figure 3. 

This observation tells us that with a spin-^ as the 
working substance, the inverse cycle of a Stirling heat en¬ 
gine may not be a refrigerator. This is totally different 
from that in classical physics. This feature is not rare, in 
fact with an ion trapped in a shallow potential, i.e., a har¬ 
monic oscillator with finite energy-levels, the inverse cycle 
of a Stirling heat engine is not a refrigerator, too. See, Fig¬ 
ures 4a-4c. However, if the trapped potential is infinitely 
high, namely, the harmonic oscillator has infinite energy 
levels, the inverse cycle is a refrigerator, numerically de¬ 
scribed as Figure 4d, which comes from the following 
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Fig. 4. Coefficient of performance £ as functions of Huj\ (in 
units of kT c ) and i- Other parameters are kTh = 3 kT c . 
The figures illustrate the heat engine composed of the working 
substance having five (a), ten (b), twenty (c), and infinite (d) 
energy levels. 


calculations: 
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3 Coupled spins 


Now we consider the coupled spins as the working sub¬ 
stance. The Hamiltonian of the working substance under 
consideration can be written as: 

# = § K + (tD + j K ^ 2 “ + 4 ) ( 15 ) 


where J is the coupling constant between the spins. J > 0 
and J < 0 correspond to the antiferromagnetic and the 
ferromagnetic case, respectively. Here we only consider the 
antiferromagnetic case, i.e., J > 0. The four eigenvalues 
and corresponding eigenstates can be easily obtained as: 


Ei = -B,: 
E 2 = - J, 

£3 = J, | 

E4 = B. 


IVt) = |00), 

hfe) = 4 (I° 1 ) - 110 )). 

^ 3 > = -h(|01) + |10)), 

^4> = 111). 




Fig. 5. Heat absorption Qh from bath at high temperature Th, 
heat release Q c to bath at low temperature T c and net work 
output as a function of coupling constant J (in units of kT c ). 
Other parameters for the cycle are kTh = 3 kT c , B\ — 3 kT c , 
B 2 / B\ = 2. 


For the coupled spins as the working substance, we take 
the coupling constant as a parameter and keep it as a 
constant during the cycle. We consider the heat engine 
cycle only. The cycle process is also shown in Figure 1 and 
the cycle direction is ABC DA. The heat exchange in each 
stage is similar to equation (6) except for a replacement 
for the internal energy and entropy as: 



B sinh B(3 + J sinh J(3 
cosh B/3 -j- cosh J(3 ’ 


S = k 


| In [2 (cosh B/3 + cosh J(3)\ + 
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As a result, the net work output per cycle can be 
written as: 


tw _ urr ^ cosh Bi(3 h + cosh J(3 h 
h ^ cosh B 2 (3 h + cosh J(3 h 


-j- kT c In 


cosh B 2 p c 3 - cosh J(3 c 
cosh Bif3 c + cosh J(3 c 



For this net work output, W > 0 and W < 0 denote the 
work is done by and on the working substance. The net 
amount of heat transfer between the system and regener¬ 
ator is AQ = Qbc + Qda- With a single spin as the working 
substance, the perfect regeneration condition AQ = 0 is 
difficult to met for arbitrary B\, £> 2 , T c , and T^. How¬ 
ever, for coupled spins, this condition can be satisfied for 
arbitrary B\, H 2 , T c and Th by adjusting the coupling con¬ 
stant J. In general, for fixed B\, B 2 and T c , T^, there are 
two coupling constants J that can enforce AQ = 0. We 
will see in the following that these two points correspond 
to two different cases. 

Next we consider the heat absorption Qh from the bath 
at high temperature T^, heat release Q c to the bath at 
low temperature T c and net work output W. As discussed 
above, when AQ ^ 0, Qh and Q c are determined by Q 1 , 
Q 2 , and AQ. The numerical results are shown in Figure 5. 
Here we calculate these thermodynamical quantities as a 
function of coupling constant J (in units of kT c ). There 
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Fig. 6. Efficiency or coefficient of performance as a function 
of coupling constant J (in units of kT c ). Other parameters in 
the cycle are: kTh = 3 kT c , B 2 /B 1 = 2. Left: B\ = kT c ; right: 
B\ = 3 kT c . The grey regimes in the figures correspond to the 
cases of neither heat engine nor refrigerator. 


are three vertical black lines in the figure, which divide 
the figure into four regimes, named 1 ~ 4. In regime 1, 
we have Qh > 0, Q c < 0, and W > 0. This is the typical 
regime for the heat engine cycle. But it is not true for 
other regimes. In regime 2, we see Qh > 0, Q c < 0 and 
W < 0. This means the system absorbs heat from bath 
Th, and releases heat to bath T c with the work done on the 
working substance. This has no application in industry. In 
regime 3, we have Qh < 0, Q c < 0 and W < 0, i.e., the 
system releases heat to both bathes and a work is done on 
the working substance. This also has no use. In regime 4, 
we obtain Qh > 0, Q c < 0 and W < 0. This is the cycle 
of refrigeration. Hence we conclude that, when increasing 
the coupling constant J, we can turn a heat engine cycle 
into a refrigeration cycle. 

The numerical results of efficiency and the perfor¬ 
mance coefficient are given in Figure 6. The two vertical 
red lines indicate AQ = 0, one of them is in the regime of 
heat engine cycle, and the other is in the regime of refrig¬ 
eration cycle. Moreover, we can see from the figure that 
when AQ = 0, both the efficiency and the coefficient of 
performance reach a peak. The peak for the coefficient of 
performance is the maximum for refrigeration cycle, how¬ 
ever, the peak for the efficiency is not always the maximal 
efficiency for the heat engine cycle. For a set of parame¬ 
ters T>i, T >2 and kT c and feT^, if the heat engine efficiency 
is lower than Carnot efficiency when J — 0, the peak is 
the maximal efficiency (Fig. 6a), otherwise it is not, see 
Figure 6b. 

As shown, the coupling constant can turn a heat engine 
into a refrigeration cycle, the physics behind this conver¬ 
sion is as follows. We plot the isothermal lines in the mag¬ 
netic field-entropy plane (B-S plane) for different coupling 
constants J. In this figure, other parameters are same as 
those in Figure 6a. There are two coupling constants cor¬ 
responding to AQ = 0 (vertical red line in Fig. 6a), which 
are marked by J\ and J 2 with J\ < J 2 . In these two cases, 
the heat released in one regenerative process is absorbed 
entirely by the other regenerative process. Comparing Fig¬ 
ures Ta and Tb, one can easily find that the monotonicity 
of the isothermal lines are different, i.e., when increasing 
the coupling constant J, the monotonicity of the isother¬ 
mal lines change. As a result, the isothermal process with 
heat absorption becomes the isothermal process with heat 




Fig. 7. The isothermal line in the magnetic field-entropy (B-S) 
plane. The parameters in the figure are same to Figure 6a. In 
Figure 6a, the two coupled constants correspond to AQ = 0 
are named as J\ and J 2 , where J\ corresponds to the maximal 
efficiency and J 2 corresponds to the maximal coefficient of per¬ 
formance. In this figure, we have set: (a) J — Ji\ (b) J = J 2 . 


release, and a heat engine cycle becomes a refrigeration cy¬ 
cle. We should note in Figure 7a that when T = T c , the 
isothermal line is a non-monotonic curve. This means in 
an isothermal process with T = T c , when we increase the 
magnetic field, the working substance absorbs heat first 
then releases heat shortly. The total result is that some 
heat is released to the cold bath. 


4 Proposed realization in a Paul trap 

The prediction can be observed in a single ion trapped 
in a shallow linear Paul trap. The shallow trap can sup¬ 
port only finite energy levels, realized by shortening the 
length of the rf electrodes in the setup in reference [31]. In 
contrary to the proposal in reference [31], we replace the 
isentropic process which is nonthermal by the isothermal 
process, this can be done by keeping the system in contact 
with the bath. By changing the length of the rf electrodes, 
we can manipulate the number of levels in the trap and 
consequently control the regime in which the inverse cycle 
of the Stirling heat engine is not a refrigerator. 

5 Conclusions 

In summary, we have considered the quantum Stirling cy¬ 
cle for a single spin or coupled spins as the working sub¬ 
stance. The quantum Stirling cycle is the quantum me¬ 
chanical generalization of the classical Stirling cycle which 
contains two quantum isothermal processes and two quan¬ 
tum isochoric processes. According to the quantum me¬ 
chanical generalization of classical thermodynamical pro¬ 
cesses, in the quantum isothermal process, the working 
substance is kept contact with a heat bath and the ex¬ 
ternal magnetic field is altered slowly so that the work¬ 
ing substance is always in the thermal equilibrium state. 
We should notice that in classical isothermal process, the 
internal energy of the working substance (the ideal gas) 
does not depend on the generalized coordinate of the sys¬ 
tem (volume V). But this is not true for quantum working 
substance. This properties leads to new result for quantum 
Stirling cycle. In the quantum isochoric process, the exter¬ 
nal magnetic field is fixed. Our results show that due to 
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the properties of quantum mechanical system, the ther¬ 
modynamical cycle exhibits many interesting properties 
which are quite different from the classical counterpart. 
For a single spin, both the heat engine cycle and refriger¬ 
ation cycle are discussed. We have found that for heat en¬ 
gine cycle, the efficiency can exceed the Carnot efficiency. 
This improvement comes from the use of the regenerator. 
The regenerator is similar to the classical. As a result, the 
reason for such a high efficiency comes from the quantum 
mechanical properties of the working substance. Moreover 
we find that the inverse cycle of the heat engine is not al¬ 
ways a refrigeration cycle. This holds true for a Stirling 
heat engine with an ion trapped in a shallow potential as 
the working substance. For the refrigeration cycle, the co¬ 
efficient of performance is always lower than its classical 
counterpart. These results indicate that quantum Stirling 
cycle is more suitable for designing a heat engine but not 
advisable for refrigerator. With two coupled spins as the 
working substance, we find that the heat engine and the 
refrigeration cycle can be altered by increasing the cou¬ 
pling constant. This alternation can be explained as the 
transformation of the monotonicity for the isothermal line 
in magnetic field-entropy plane. 


This work is supported by NSF of China under Grant 
Nos. 11105064 and 11175032. 


Appendix A: The proof for the positive 
work condition 


Here we give the detailed proof for the positive work 
condition. We first define a function as: 


cosh(x /2kTh)] kTh 
cosh (x / 2 kT c ) ] kTc 


(A.l) 


Taking derivative with respect to x, we have: 


, 1 [cosh(x/2 kT h )] kTh 1 . (kT c - kT h )x 

j (x) = - —— , sinh 


2 [cosh(x/2 kT c )\ 


kT c -\-l 


2 kT c kTh 


(A.2) 


Here we assume that T c < Th , as a result, for V x > 0, 
we have f'(x ) < 0. This means that f(x) is a decreasing 
function. Consequently, if B\ < B 2 we obtain f(B 1 ) > 
f(B 2 ). This inequality leads to W >0. 


The first step A —> B is an isothermal process, in 
which dQ = TdS. As a result 


>B 


A 


dQ 

T 


f dS = S(B) - S(A) 

Ja 

/ 1 u PhBi kfdhBi PhB\ 

k In cosh-tanh 


7 1 u , kf3 h B 2 

— k In cosh-1-tanh 


2 

PhB 2 


(B.2) 


The second step B —> C is an isochoric process, in which 
we have dQ = E n dp n - We can obtain: 


*c 


B 


dQ 

T 


■c 




(B.3) 


For a single spin, the two eigenvalues and the corre¬ 


sponding probabilities are A ,2 = 1 and pi, 2 


kT sech^f^, respectively. Based on these results, we 
have: 


>c 


B 


dQ 

T 


4 k 


* T r 


T h 


— k In cosh 


sec h 2 A. 

,‘ 2kT dT 

(3 c Bi kflcB! j3 c B\ 

-tanh- 


- k In cosh Ai + tanh Ai. (B.4) 


The following two steps can be calculated similarly as: 


,D dQ 

c T 


ii : kf3 c B2 dcB‘> 

k In cosh-tanh- 


k In cosh Ai + MTl tanh BA ( B .5) 


• A 


D 


dQ 

T 


1 1 U PhB 2 k(3 h B 2 
k In cosh-tanh 


2 

PhB 2 


k In cosh + k ^ c ^ 2 tanh . (B.6) 


Combining equations (B.2)-(B.6), we have f ^ = 0, i.e., 
the quantum Stirling cycle with single spin as working 
substance is reversible. 


Appendix B: Proof of the the reversibility 
of the quantum Stirling cycle with single 
spin as working substance 

According to the second law of thermodynamics, a 
reversible cycle should satisfy 



in a whole cycle. We consider one direction for the cy¬ 
cle, for example, heat engine cycle ABC DA in Figure 1. 
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